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NUMERICAL STUDY OF THE STABILITY OF THE PEREGRINE 

BREATHER 

CHRISTIAN KLEIN AND MARIANA HARAGUS 


Abstract. The Peregrine breather is widely discussed as a model for rogue waves in deep 
water. We present here a detailed numerical study of perturbations of the Peregrine breather 
as a solution to the nonlinear Schrodinger (NLS) equations. We first address the modulational 
instability of the constant modulus solution to NLS. Then we study numerically localized and 
nonlocalized perturbations of the Peregrine breather in the linear and fully nonlinear setting. It 
is shown that the solution is unstable against all considered perturbations. 


1. Introduction 


The importance of rogue waves needs hardly be stressed, and the attempts in the mathematical 
modelling of these phenomena have been intensified in recent years. Since the governing equations 
for hydrodynamics are similar to the ones from nonlinear optics, also optical rogue waves have 
been experimentally and analytically studied recently, see [15] for a review. In deep water and 
nonlinear optics the focusing nonlinear Schrodinger (NLS) equation 

( 1 ) iut+Uxx + ‘2.\u\^u = Q. 


for a complex u{x, t) describing the modulation of a carrier wave is an important model equation. 
This is because its solutions show nonlinear focusing and a modulation instability, i.e., self-induced 
amplitude modulation of a continuous wave propagating in a nonlinear medium, with subsequent 
generation of localized structures. For a discussion of NLS in the context of rogue waves see for 
instance [20] and references therein. It is fair to say that no complete understanding of rogue 
waves has yet been reached. For complementary approaches see for instance [5, 6]. 

Remarkably the NLS equation has solutions called breathers that show the features of the 
popular definition of rogue waves as ‘waves that appear from nowhere and disappear without a 
trace’, see for instance [1]. The most popular nonperiodic breather solution in this context is the 
Peregrine breather [21] 


( 2 ) 


Uper 


/ _ 4(1 -h 4it) \ 2it 

V 1 + 4x2 + 16t2 J 


which is an exact solution of the focusing NLS equation. It has the property to have constant 
modulus for ]<] —>■ oo, thus appearing from nowhere and disappearing after some time, and for 
|x| —>■ oo. The maximum of |u| is reached for x = t = 0 and is three times the asymptotic value. 
Exact NLS solutions describing multiple breathers are given in [12], for an alternative approach to 
generate breather solutions see [16]. It is reported that the Peregrine solution has been recently 
observed experimentally in hydrodynamics [8, 9], in plasma physics [2] and in nonlinear optics 

[17]. 

For a solution to be of practical relevance in the context of rogue waves, it must be sufficiently 
robust to be observable in a realistic environment where the initial conditions are never exactly 
given by the solution (2), but are at best random perturbations thereof. This means one has to 
study perturbations of Peregrine initial data and to investigate whether the resulting perturbed 
solution to the NLS equation is in some sense close to (2). This is a difficult task mathematically 
even in a linearized setting, i.e., for small amplitude perturbations, since the background solution 
is not constant in time. The question of nonlinear stability is even more out of reach. Thus 
with the currently available techniques, one has to address this question numerically. For periodic 
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breathers, this has been recently done in [7] . An experimental study of the stability of the Peregrine 
breather appearing in hydrodynamical experiments against wind was described in [10]. 

It is the goal of the present paper to present a detailed numerical study of the stability of the 
Peregrine solution. We consider first small perturbations, e.g., the evolution of initial data of 
the form u{x,0) = uper{x,0) + v{x) with ||u(a;)||oo 1 according to the NLS equation (1). Of 

special interest are localized perturbations v(x) in the Schwartz space of smooth rapidly decreasing 
functions S(M). The second kind of perturbation to be studied is of the form uq(x) = auper{x, to) 
with (7^1. This means we will study both rapidly decreasing and non decreasing perturbations 
of the Peregrine breather. 

The paper is organized as follows: in section 2 we collect some mathematical preliminaries 
needed for the following sections. In section 3 we study the absolute instability of the asymptotic 
state of the Peregrine solution, i.e., of a solution with constant modulus. In section 4, we present 
a numerical study of the linearized (on the Peregrine background) NLS equation for localized and 
nonlocalized perturbations. In section 5, localized and nonlocalized perturbations of the Peregrine 
breather are studied for the full NLS equation to explore the nonlinear stability of the breather. 
We add some concluding remarks in section 6. 


2. Preliminaries 


In this section we collect some mathematical information on the NLS equation and linearizations 
on the background of the Peregrine breather needed for the ensuing stability study. We also briefly 
review the numerical techniques to be applied in the following. 


2.1. NLS and linearized equation. It was shown by Zakharov and Shabat [23] that the NLS 
equation is completely integrable, i.e., it has an infinite number of conserved quantities the most 
prominent of which being the mass, the square of the norm of u and the energy. Since we want 
to discuss in this note solutions not decaying to 0 at infinity but with lim| 2 ,|_).oo [m] = k > 0 (i.e., 
solutions which are not in L^), we consider instead the quantity 



which is a combination of mass and energy. The complete integrability of the NLS equation is the 
reason why many explicit solutions as solitons or breathers are known. 

The NLS equation has also some conformal symmetry which implies that with u{x,t) also 

(4) u‘^(x,t) = au(ax,a^t), cr G K, 


is a solution of NLS. This mean that asymptotically nonvanishing NLS solutions could be always 
transformed by a suitable choice of a in (4) to have modulus 1 at infinity. 

We first study linearizations of the NLS equation (1) of the form u = uper + v which leads in 
lowest order of v to the equation 

(5) iVt+Vxx+^\uper\'^V + 2u^p^.^V = Q. 

It is straight forward to check that this equation has the conserved quantity (this is a consequence 
of the conservation of the norm for solutions to the NLS equation) 


(6) M = J [uperV + UperV)dx. 

Note that there are different ways to linearize the NLS equation (1). An alternative form is to 
write u = uper(l + v) which implies 

(7) ivt + Vxx + 2{\nUper)xVx + ‘i\uper\‘^^V = 0, 


for NLS (1) in lowest order of v. The disadvantage of equation (7) from a numerical point of view 
is that it is singular at the zeros of the Peregrine solution for t = 0. Even for finite small t, the 
fact that one is close to a singularity is numerically problematic. Therefore we study in this note 
for linearizations numerically only equation (5). However, equation (7) is convenient for analytical 
approaches, see the following section. 
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2.2. Numerical methods. We use in this paper spectral methods in a;, i.e., methods which show 
for a smooth solution an exponential decrease of the numerical error with the number of collocation 
points. For the time integration, we use fourth order methods. The accuracy of the numerical 
computation is controlled via conserved quantities. Because of unavoidable numerical errors, these 
are in an actual computation not exactly conserved, but can be used as a valid indicator of the 
reached precision: as shown for instance in [19], they tend to overestimate the numerical accuracy 
by 1-2 orders of magnitude. 

For rapidly decreasing v{x, •) in (5), we use a Fourier spectral method for the spatial dependence 
and an exponential time differencing (ETD) scheme due to Cox and Matthews [11] for the time 
integration. The latter scheme uses an exact integration of the linearized equation in Fourier space. 
The modulus of the Fourier coefficients in the spatial coordinate indicates the spatial resolution 
since it is known that they decrease exponentially for analytic functions. In all studied examples 
in this paper, they decrease to machine precision (here « 10“^®) during the whole computation 
which indicates that the solution is always fully resolved spatially. 

For nonlocalized perturbations of the Peregrine solution (2) in the full NTS equation (1), Fourier 
methods are not appropriate. The algebraic decrease of the solution for |a:| —> oo makes it difficult 
to periodically continue the solution as an analytic function within the finite numerical precision 
as can be done with Schwartz functions. To obtain spectral convergence also in this context, a 
novel numerical approach has been presented in [4]: the real line is divided into several domains 
each of which is mapped to [—1,1] with a Mobius transformation. Near infinity we use 1/x as 
a local coordinate which allows to cover the whole real line with this approach. Note that we 
use in contrast to [4] four domains here among which just one compactified one (with 1/a: as a 
coordinate). This is due to the symmetric (in x) problem studied here in contrast to [4]. On each 
domain we use polynomial interpolation to obtain on each interval a spectral method. Imposing 
at the domain boundaries the conditions that the solution be there implies with the analytic 
properties of the NTS equation (which is a second order PDF in x) that we obtain a spectral 
method on the whole real line. For the time integration, we use again a fourth order method. 

For details the reader is referred to [19] and [4] and references given therein. Note that the 
numerical approaches are completely independent for the localized perturbations in the linear case 
and the nonlinear case. Therefore we compare the results in these cases showing that they are 
qualitatively similar until the nonlinearity dominates in the latter case. This gives further evidence 
for the validity of the presented results. 

3. Absolute instability of the asymptotic state 

A problem in the stability study of the Peregrine solution is the fact that it depends both 
on X and t. Therefore, for the analytical study in this section, we focus on the stability of the 
asymptotic state for x —oo when uper —t Instabilities of asymptotic states typically induce 

instabilities of the corresponding exact solutions, however we do not attempt to give a rigorous 
proof of this fact for the Peregrine solution. 

The asymptotic state of the Peregrine solution uper has constant modulus, so that for our 
analysis it is more convenient to use the linearized equation (7). For complex perturbations of the 
form V = a + ip we obtain the system 

(8) “b Pxx — 0, Pt ^xx 4oi = 0, 

which has constant coefficients. Our purpose is to study the spectral properties of the matrix- 
operator defined through 

l3(o;, /3) — ( Pxx^ ^xx 4“ do), 

and acting in (M) x (M). This function space corresponds to a choice of localized perturbations, 
but the instability results below also hold for the larger class of nonlocalized perturbations. We 
compute the spectrum, which coincides with the essential spectrum, of the operator £, and then 
the absolute spectrum. Absolute spectra, which are always located to the left of essential spectra 
in the complex plane, allow to distinguish between absolute instabilities, when perturbations grow 
in time at each point of the domain, and convective instabilities, when though overall norms grow 
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in time, perturbations locally decay: growing perturbations are convected away towards infinity 
(e.g., see [22] and the references therein). 

Since the operator C has constant coefficients, these spectra can be determined from the dis¬ 
persion relation 

(9) V{\v) = X^ + + 

obtained by looking for solutions of (8) of the form {a, /3){x, t) = b), for complex numbers 

A, V and a, b. First, the essential spectrum, which coincides with the full spectrum of the operator 
C in this case, is the set 

CT = {A e C ; I?(A, ik) = 0, for some k G M}. 

Then 

cr = {A G C ; A^ = fc^(4 — /c^), for some k G M}, 

implying ct = U [—2, 2]. In particular, the set cr contains values in the open right hand complex 
plane showing that the asymptotic state of the Peregrine solution is (essentially) spectrally stable. 

Next, the absolute spectrum is found by looking at the relative location in the complex plane 
of the solutions v G C oi the dispersion relation = 0 for A G C [22]. More precisely, for 

complex numbers A with large real part, the dispersion relation (9) has two roots i4+(A), 1^2 W 
with positive real parts, Rei/(''(A) > Rer'^(A) > 0, and two roots i^]~(A), r'^(A) with negative real 
parts, Rei/]"(A) < Rei/^(A) < 0. Then in this case, the absolute spectrum of the operator £ is 
simply defined by 

CTabs = {A G C ; Rei/^(A) = Rei/^(A)}. 

The essential instability found above is absolute if there exists A G cTafcs with Re A > 0. Since 
the dispersion relation is invariant under the change 1 / 1 —> — 14 , we have that vJ{X) = (A), for 

j = 1,2, implying that Rez/^(A) = Re^'^(A) = 0 when A G (Jabs- As a consequence, the absolute 
spectrum coincides with the essential spectrum here, aabs = cr = U [—2,2], showing that the 
instability of the asymptotic state of the Peregrine solution is absolute. 

4. Numerical simulations for the linearized equation 

In this section, we study perturbations of the Peregrine breather in the linearized regime. This 
means we study initial value problems for the linear equation (5) both for rapidly decreasing initial 
data and for data proportional to the Peregrine solution. We find in both cases that the solution 
grows strongly in time which confirms the absolute instability discussed in the previous section. 
Because of the time dependent background leading to a nonautonomous equation (5), the growth 
in time is not exponential. 

4.1. Rapidly decreasing perturbations. For the numerical experiments, we concentrate on 
equation (5) and the initial data v = O.le”"^ . Since these data are rapidly decreasing, we use as 
discussed in section 2 a Fourier spectral method for the spatial dependence. The computation is 
carried out with N = 2^^ Fourier modes and Nt = lO'^ time steps. The modulus of the Fourier 
coefficients decreases to machine precision during the whole computation which indicates that the 
solution is fully resolved spatially. The relative change of the conserved quantity (6) during the 
computation is of the order of 10“^° which indicates that the solution shown in Fig. 1 should be 
accurate at least to the order of 10“®. The maximum of the solution is visibly quickly growing 
and becomes larger than the modulus of the Peregrine solution itself which indicates that the 
linearized regime is left, in other words that the solution is unstable. Moreover the initial peak is 
also being dispersed to ±00 as can be seen in Fig. I. 

4.2. Nonlocalized perturbations. We consider nonlocalized data proportional to the breather. 
We use the same numerical approach as for the fully nonlinear equation and the same parameters 
as detailed there. For the initial data v(x, 0) = 0.1uper(a:, 0), we get the solution shown in Fig. 2. 
As for the localized perturbations in Fig. 1, the solution is clearly unstable and quickly leaves the 
linearized regime. 
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Figure 1. Numerical solution to the linearized NLS equation (5) for the initial 
data V = 0.1e“^^ 



Figure 2. Numerical solution to the linearized NLS equation (5) for the initial 
data V = Q.luper 


5. Nonlinear Stability 

In this section, we will study perturbations of the Peregrine breather for the full NLS equation 
(1). We again consider localized and nonlocalized perturbations to the Peregrine breather as in the 
linearized case and find that the breather is unstable against all perturbations. This also applies if 
the solution is perturbed at an earlier time: the resulting NLS solution for these perturbed initial 
data will not be close to the Peregrine breather at t = 0 or later times. 
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The numerical approach uses four domains for x: domain I with x S [—20, —5], domain II with 
X G [—5,5], domain III with x G [5, 20]; the compactified domain IV is defined via —1/20 < 1/x < 
1/20. Thus we solve the NTS equation on the whole real line. In the respective domains, we use 
Ni = Niii = 200, Nil = 300 and Njv = 100 Chebyshev polynomials. We always take the time 
step At = 10“^. The numerical accuracy is controlled via the conserved quantity £ in (3) in the 
form of the relative quantity 


( 10 ) 


Ae = 1 — 


m 

sitoY 


(because of numerical errors the computed £ will be time dependent). 


5.1. Rapidly decreasing perturbations. We first study as in the previous section initial data 
of the form u{x, 0) = uper{x, 0) + 0.1 exp(—a;^). The resulting solution can be seen in Fig. 3 where 
the initial hump appears to be dispersed away to infinity, but in a significantly different way than 
the Peregrine breather. 



Figure 3. Solution to the NTS equation for the initial data uper{x,0) + 
0.1exp(—x^) in dependence of time. 

This is even more obvious in the left figure of Fig. 4 where the solution at the last computed time 
is shown together with the Peregrine solution (2) at the same time. Clearly the Peregrine breather 
is unstable against this type of perturbations, and the perturbed solution does not stay close to 
the exact solution. Note that the solution is even at the last computed time well resolved spatially 
as can be seen in the right figure of Fig. 4, where the Chebyshev coefficients at the last computed 
time are shown. Since the initial data are symmetric with respect to the transformation x —>■ —x 
and since the Schrodinger equation preserves parity, the Chebyshev coefficients in zone I and III 
are identical (therefore the coefficients in zone III are not shown), and half of the coefficients in 
zone II and IV vanish with numerical precision. The coefficients in zone IV still decrease to 10“^. 
As discussed in [4], the slow decrease of the Chebyshev coefhcients is due to the oscillations with 
a slowly decreasing amplitude for jxj —>■ oo. The conserved quantity (10) is Ap ~ 8.8 * 10“^. 
This implies together with the resolution in coefficient space shown in Fig. 4 that the solution is 
computed to at least plotting accuracy. 

In Fig. 5 we show the difference between the numerical solution u in Fig. 3 and the Peregrine 
solution. This corresponds to the function v in the linearized case (5), and Fig. 5 thus has to be 
compared to Fig. I. It is interesting to see that both figures are qualitatively similar, but that the 
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Figure 4. Solution to the NFS equation shown in Fig. 3 at t = 1 in blue and 
the Peregrine breather (2) for t = 1 in green on the left, and the corresponding 
Chebyshev coefhcients for the numerical solution on the right (in blue for domain 
I, green for domain II and red for domain III). 

nonlinearity has the effect to delimit the growing of |u — uper\ for later times and to disperse the 
solution more efficiently. 



Figure 5. Difference between the solution to the NFS equation shown in Fig. 3 
and the Peregrine breather (2). 

It is an interesting question whether a perturbation of the breather at an earlier time still 
allows to observe a Peregrine type behavior at t ~ 0 at least in a qualitative sense. To answer 
this question, we consider as initial data u(a;, —I) = uperix, —1) + 0.1exp(—a;^) in Fig. 6. It can 
be recognized that the solution grows to a similar maximal value as the Peregrine solution, but 
reaches it at an earlier time. The focusing effect of the focusing NFS equation is clearly visible, 
but it does not force the solution towards the Peregrine breather. 

In Fig. 7, it can be seen that the solution at the last recorded time is clearly different from 
the Peregrine solution 2 which again confirms that the solution is unstable. The Chebyshev 
coefficients at the last recorded time can be seen in Fig. 7 on the right showing that the solution 
is well resolved. 
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Figure 6. Solution to the NLS equation for the initial data M(a;, —1) = 
uper{x, —1) + 0.1exp(—in dependence of time. 




X 


Figure 7. Solution to the NLS equation shown in Fig. 6 at t = 0 in blue and 
the Peregrine breather (2) for t = 0 in green on the left, and the corresponding 
Chebyshev coefficients for the numerical solution on the right (in blue for domain 
I, green for domain II and red for domain III). 


5.2. Nonlocalized perturbations. In this subsection we study nonlocalized perturbations of 
the breather of the form it(x,to) = o’Mper(a^, to)i where cr ~ 1. The parameters for the numerical 
computation are as in the previous subsection. 

In Fig. 8 we show the solution for the initial data u{x,0) = l.luper(a;, 0) in dependence of 
time. It can be seen that the solution keeps growing despite being initially already above the 
global maximum of the Peregrine solution before decreasing, but in contrast to the latter not 
monotonically. Instead it grows again before being eventually dispersed. Note that the quantity 
£ in (3) is positive in this case. 

It is obvious, see also the left figure in Fig. 9, that the solution at the final recorded time is not 
close to the Peregrine solution at the same time. Clearly the Peregrine solution is also nonlinearly 
unstable against unbounded perturbations proportional to the solution itself. The right figure in 
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Figure 8. Solution to the NLS equation for the initial data u{x,0) = 
l.luper{x,0) in dependence of time. 


Fig. 9 shows that the solution is again well resolved spatially even at the final time of computation. 
For the quantity (10) we have Ap < 8.15 * 10“"^ during the whole computation. 




Figure 9. Solution to the NLS equation shown in Fig. 8 at t = 1 in blue and 
the Peregrine breather (2) for t = 1 in green on the left, and the corresponding 
Chebyshev coefficients for the numerical solution on the right (in blue for domain 
I, green for domain II and red for domain III). 

In Fig. 10 we show the difference of the numerical solution in Fig. 8 and the Peregrine breather 
in dependence of t. This can be compared to the solution of the linearized equation (5) in Fig. 2. It 
can be seen that the solution to the linearized equation (5) grows continually, whereas the solution 
to the full NLS equation reaches a maximum and decreases then. Note however the qualitative 
similarity at earlier times. 

A similar perturbation as in Fig. 8 but of positive £ is considered in Fig. 11 where the initial 
data are u{x,0) = 0.9uper(a;, 0). Here the solution never reaches the maximum of the Peregrine 
solution. Instead it decreases in modulus with time, but also not as the latter solution, but by 
radiating the mass towards infinity in the form of dispersive oscillations. 
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Figure 10. Difference of the solution to the NFS equation shown in Fig. 8 and 
the Peregrine breather (2). 



Figure 11. Solution to the NLS equation for the initial data u(x,0) = 
0.9uper(x,0) in dependence of time. 


Fig. 12 shows that the solution is again in no sense close to the Peregrine solution (2). The 
Chebyshev coefficients in the same figure decrease again to the order of 10“^. During the whole 
computation we have Ap < 4.43 * 10“®. 

In Fig. 13 we show the difference of the solution in Fig. 11 and the Peregrine breather which 
can be compared to the solution of the linearized equation (5) in Fig. 2. Note that the solution of 
the linearized equation (5) is identical for the cases corresponding to Fig. 8 and Fig. 11. Clearly 
this is not the case for the solutions to the full NLS equation, but there is good agreement for 
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Figure 12. Solution to the NLS equation shown in Fig. 11 at t = 1 in blue and 
the Peregrine breather (2) for t = 1 in green on the left, and the corresponding 
Chebyshev coefficients for the numerical solution on the right (in blue for domain 
I, green for domain II and red for domain III). 

early times. For larger times, the difference in Fig. 13 grows considerably more slowly than the 
solution in Fig. 2. 



Figure 13. Difference of the solution to the NLS equation in Fig. 11 and the 
Peregrine breather. 

As in the previous subsection, we also address the question whether perturbations of the above 
kind at an earlier time still allow to qualitatively observe characteristic features of the Peregrine 
breather (2) at later times. In Fig. 14 we consider the solution for initial data of the form 
u{x, —1) = l.luperix, —1). In this case £ is again positive, and a higher maximum is reached at 
an earlier time than for the Peregrine breather, though the solution is qualitatively of a similar 
form. However, the solution is dispersed away to infinity in a completely different manner than 
for the Peregrine breather. 

In Fig. 15 it is, however, obvious that the solution at the final recorded time is not close to the 
Peregrine breather. The Chebyshev coefficients in the same figure decrease to 10“®. During the 
whole computation we have < 2.5 * 10“^. 
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Figure 14. Solution to the NLS equation for the initial data u(a;, —1) = 
l.luper(x, —1) in dependence of time. 




Figure 15. Solution to the NLS equation shown in Fig. 14 at t = 1 in blue and 
the Peregrine breather (2) for t = 1 in green on the left, and the corresponding 
Chebyshev coefficients for the numerical solution on the right (in blue for domain 
1, green for domain 11 and red for domain Ill). 

For the initial data u(ai, —1) = 0.9uper{x, —1), i.e., a perturbation of negative £ in (3), the 
solution is shown in Fig. 16. It can be seen that the solution does not reach the maximum of the 
Peregrine breather. 

Again the solution is not close to the Peregrine breather at the final recorded time as can be 
seen in Fig. 17. The solution is well resolved in the space of Chebyshev coefficients as can be seen 
in the same figure. During the whole computation we have < 9.05 * 10“^. 

6. CONGLUSION 

The results of the previous section indicate that solutions to the focusing NLS equation are not 
stable if they do not vanish for |a;| —?► oo, and that this instability is absolute. This appears to 
be the case whether the perturbation is localized or not. On the other hand the focusing effect 
of the NLS equation is of course always present. In fact it seems that the Peregrine breather is 
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Figure 16. Solution to the NLS equation for the initial data u(a;, — 1) = 
0.9uper{x, —1) in dependence of time. 




Figure 17. Solution to the NLS equation shown in Fig. 16 at t = 1 in blue and 
the Peregrine breather (2) for t = 1 in green on the left, and the corresponding 
Chebyshev coefficients for the numerical solution on the right (in blue for domain 
1, green for domain 11 and red for domain Ill). 


not stable, but appears in an asymptotic sense explained below locally in strongly focused peaks. 
Note that the behavior shown above is actually observed in experiments in nonlinear optics for 
initial data close to the Peregrine breather [18]. 

If one considers initial data with a support on scales of order 1/e and studies the corresponding 
NLS solution on time scales of order 1/e, this can be conveniently done by the coordinate change 
X ^ ex, t ^ et which leads for (1) to 

( 11 ) ieut + e^Uxx +‘^\u\^u = 0- 

It is well known that solutions to the NLS equation in the semiclassical limit e ^ 1 for analytic 
initial data in with a single maximum become strongly peaked at t ^ tc where tc is the critical 
time where the solution to the semiclassical system (the system following from (11) by formally 
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letting e —> 0) for the same initial data develops a cusp. An example for Gaussian initial data is 
shown in Fig. 18. 



Figure 18. Solution to the NFS equation (11) with e = 0.1 for the initial data 
uq = exp(—cc^). 

In [13, 14] it was conjectured that the local (in the vicinity of the cusp) behavior of the NFS 
solution for t tc is asymptotically for e —>■ 0 given by the tritronquee solution to the Painleve 
I equation. A partial proof of this conjecture was given in [3]. In [3] it was also shown that the 
Peregrine breather appears in the asymptotic (local) description of the peak of the NFS solution 
in the semiclassical limit and is in this sense universal. 
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